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ABSTRACT 

This  is  the  second  part  of  a  study  which  deals  with  the  problem  of 

passive  time  delay  estimation.  The  focus  here  is  on  systems  employing 

wideband  signals  and/or  arrays  of  very  widely  separated  receivers.  A 

modified  (improved)  version  of  the  Ziv-Zakai  lower  bound  (ZZLB)  is  used  to 

analyze  the  effect  of  additive  noise  and  signal  ambiguities  on  the  attainable 

mean  square  estimation  errors.  When  the  lower  bound  is  plotted  as  a  function 

of  signal-to-noise  ratio  (SNR)  one  observes  two  distinct  threshold  phenomena 

dividing  the  SNR  domain  into  three  disjointed  segments:  at  high  SNR  the 

lower  bound  coincides  with  the  Cramer-Rao  lower  bound  (CRLB).  This  is  the 

ambiguity-free  mode  of  operation  where  differential  delay  estimation  is 

subject  only  to  local  errors.  At  moderate  SNR  (between  the  two  thresholds), 
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the  lower  bound  exceeds  the  CRLB  by  a  factor  of  12(a)0/W)  where  <dq  and 
W  are,  respectively,  the  center  frequency  and  signal  bandwidth.  In  this 
region  the  ambiguities  in  the  received  signal  phases  cannot  be  resolved, 
however  a  useful  estimate  of  the  differential  delay  can  still  be  obtained 
using  the  received  signal  envelopes.  At  low  SNR,  the  lower  bound  approaches 
a  constant  level  depending  only  on  the  variance  of  the  a-priori  search  domain 
of  the  unknown  delay  parameter.  In  this  region  signal  observations  are 
subject  to  envelope  ambiguities  as  well,  thus  essentially  useless  for  the 
delay  estimation. 


I.  INTRODUCTION 
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A.  Ambiguity  Phenomena  in  Time  Delay  Estimation 
Estimation  of  the  time  difference  of  arrival  of  a  noise-like  random 
signal  observed  at  two  or  more  spatially  separated  receivers  is  a  problem  of 
consiaerable  practical  interest  in  many  disciplines  such  as  underwater 
acoustics,  geophysics  and  radio  astronomy,  to  mention  a  few.  Consequently, 
numerous  procedures  have  been  proposed  for  passive  time-delay  estimation  (e.g. 

In  most  of  the  systems  which  have  been  analyzed  the  Cramer-Rao 
inequality  has  been  used  to  set  a  lower  bound  on  the  attainable  mean  square 
estimation  errors.  Its  use  was  justified  by  invoking  a  well  known  theorem  in 
statistics  asserting  that  the  maximum  likelihood  (ML)  estimator  is 
asymptotically  unbiased  and  that  its  error  variance  approaches  the  Cramer-Rao 
lower  bound  (CRLB)  arbitrarily  close  for  sufficiently  long  observation  times. 
There  remains  the  question:  How  long  is  "long  enough"?  Clearly,  the 
observation  time  T  must  be  large  compared  with  the  correlation  time  (inverse 
bandwidth)  of  signal  and  noise  (WT/2*  >>  1),  a  condition  which  presents  very 
little  difficulty  in  practice.  However,  if  one  examines  the  problem  more 
closely,  one  finds  a  second  condition:  the  ML  estimator  (which  approaches  the 
CRLB  asymptotically)  must  not  be  subject  to  ambiguities. 

Perhaps  the  most  common  setting  in  which  this  difficulty  occurs  is  in 
aelay  estimation  using  very  narrowband  signals.  Consider  the  extreme  case  of 
observations  at  only  two  receivers  so  that  only  one  differential  delay  can  be 
estimated.  The  ML  estimate  of  that  delay  cross  correlates  the  received 
signals,  averages  for  time  T,  and  obtaines  the  desired  delay  from  the  peak  of 
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the  cross  correlation  function.  In  the  narrowband  case,  the  differential 
aelay  causes  essentially  to  a  phase  shift  between  the  received  signals  and 
thus  generating  a  formidable  ambiguity  problem.  This  is  illustrated  in  Figure 
1.  The  cross  correlation  output  peaks  at  ax,  the  true  differential  delay,  but 
is  quasi-periodic  with  a  period  of  2Wu>0  where  tu0  is  the  center  frequency 
of  tne  signal.  To  come  close  to  the  CRLB  one  must  be  able  to  distinguish 
unambiguously  between  adjacent  peaks  of  the  correlation  function.  If  the 
signal  bandwidth  is  only  a  small  fraction  of  its  center  frequency  (i.e., 

W/o)0  «  1),  adjacent  peaks  have  very  nearly  equal  height  and  identification 
of  the  largest  one  will  require  either  very  large  SNR  or  exceedingly  long 
observation  times.  In  many  important  practical  situations,  therefore,  the 
attainable  m.s.e.  may  be  very  drastically  inferior  to  that  predicted  by  the 
CRLB. 

In  [8 j ,  a  new  lower  bound,  based  on  a  modified  (improved)  version  of  the 
Ziv-Zakai  lower  bound  (ZZLB),  is  developed  to  analyze  the  attainable  mean 
square  error  (m.s.e.)  in  delay  estimation  schemes.  The  resulting  lower  bound 
is  then  applied  to  investigate  the  effect  of  ambiguity  on  delay  estimation 
using  narrowband  signals  [9].  When  the  lower  bound  is  plotted  as  a  function 
of  SNR,  one  observes  a  strong  threshold  phenomenon.  Above  a  critical  SNR  the 
lower  bouna  approaches  the  CRLB.  In  this  region  the  indicated  phase 
ambiguities  in  the  differential  delay  observations  can  essentially  be 
resolved.  Below  the  threshold,  the  lower  bound  exceeds  the  CRLB  by  a  large 
factor.  The  point  at  which  the  threshold  occurs  as  well  as  its  magnitude  are 
factors  of  obvious  practical  interest  since  system  performance  predictions 
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generally  are  based  on  the  small  error  assumption  (i.e.,  CRL8)  and  may, 
therefore,  be  extrapolated  incorrectly  below  the  threshold  yielding  extremely 
optimistic  performance  predictions. 

In  the  second  part  of  the  study  we  broaden  the  range  of  investigation  to 
wideband  signals  which,  for  the  delay  estimation  problem,  will  be  divided  into 
two  categories:  baseband  signals  and  bandpass  signals. 

In  the  basebana  case,  a  typical  cross  correlation  response  is 
illustrated  in  Figure  2  where  W  denotes  signal  bandwidth  and  [-D/2,  D/2] 
denotes  the  a-priori  search  domain  of  the  unknown  delay  parameter.  One 
immediately  observes  a  similar,  though  not  as  critical,  ambiguity  phenomenon 
whenever  the  a-priori  search  domain  contains  at  least  several  peaks  of  the 
cross  correlation  function  (i.e.,  WD/2ir  >>  1). 

In  the  bandpass  case  a  typical  cross  correlation  response  is  illustrated 
in  Figure  3.  Here  we  observe  two  types  of  ambiguities:  the  more  critical 
ambiguity  problem  results  from  the  highly  oscillatory  nature  of  the  phase  of 
the  observed  cross  correlation  function.  A  secondary  ambiguity  phenomenon 
results  from  the  oscillatory  nature  of  the  envelope  of  the  cross  correlation 
fu notion. 

The  purpose  of  this  paper  is  to  analyze  the  threshold  effects  and  the 
attainable  m.s.e.  for  these  two  classes  of  signals. 

B.  Problem  Formulation  and  Assumptions 

The  basic  system  of  interest  here  consists  of  a  stationary  source 
radiating  a  noise-like  signal  towards  two  spatially  separated  receivers.  Each 
receiver  also  receives  an  additive  noise  component  so  that  the  actual 
waveforms  observed  at  the  receiver  outputs  are  given  by 
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r^t)  =  s  ( t )  +  n^(t) 

-  T/2  <  t  <  T/2  (1) 

r2(t)  =  s(t  -  at)  +  n2(t) 

We  shall  assume  that  s(t),  n^(t)  and  n2(t)  are  sample  functions  from 
uncorrelatea  zero-mean  Gaussian  random  processes  with  spectral  densities  S(w), 
N^(u)  and  N2(u),  respectively.  Since  we  are  primarily  interested  in  the 
ambiguity  problem  in  time  delay  estimation,  implying  the  use  of  widely 
separated  receivers,  the  assumption  of  noise  incoherence  from  receiver  to 
receiver  is  likely  to  be  satisfied. 

We  shall  further  assume  that  at,  the  receiver-to-receiver  delay,  is 
uniformly  distributed  in  the  interval 

-D/2  <  at  <  D/2  (2) 

This  a-priori  domain  may  come  about,  perhaps  from  the  known  receiver 
separation  and  the  known  velocity  of  propagation  in  the  medium. 

Finally,  we  shall  assume  that  the  observation  time  T  is  large  compared 
with  the  correlation  time  (inverse  bandwidth)  of  signal  and  noise,  i.e. 

WT/2ir  »1.  This  condition  is  very  generally  satisfied  in  problems  of 
practical  interest  here. 

The  problem  may  now  be  stated  as  follows:  given  the  data  at  the 
receiver  outputs  (i.e.,  r^(t)  i  *  1,2),  characterize  the  minimum  m.s.e. 
estimate  of  the  delay  parameter.  Our  approach  is  to  set  a  lower  bound  on  the 
attainable  m.s.e.  using  the  modified  ZZLB.  This  approach  is  completely 
independent  from  the  actual  estimation  method.  However,  in  [10]  it  is  shown 
that  for  sufficiently  large  WT  products  the  cross  correlator  performance  comes 
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close  to  the  lower  bound  below  as  well  as  above  the  threshold.  These  results 
establish  the  modified  ZZLB  as  an  extremely  tight  lower  bound  for  problems  of 
this  type  while  demonstrating  that  the  cross  correlator  is  a  nearly  optimal 
instrumentation  in  both  the  small  and  large  error  regimes. 

C.  The  Modified  Ziv-Zakai  Lower  Bound 

The  modified  ZZLB  on  the  m.s.e.  of  any  estimate  ax  of  at  is  given  by 
(L8L  19J) 

_>-jj  f  xG[(D-x)Pe(x)]  dx  (3) 

0 

where  G[  ]  is  a  non-increasing  function  of  x  obtained  by  filling  the  valleys 
in  the  function  (D-x)Pe(x)  (see  Figure  B-l  in  Appendix  B).  Pe(x)  is  the 
minimum  probability  of  error  (achievable  by  the  likelihood  ratio  test)  for 
deciding  whether  the  true  value  of  the  parameter  is  axQ  or  ax^  where  Aij 
-  axq  =  x.  In  general,  a  closed  analytical  form  for  Pe(x)  cannot  be  found. 
However,  in  ([9],  Appendix  A)  it  has  been  shown  that  for  WT/2n  >>  1,  Pe(x)  is 
very  closely  approximated  by 

Pe(x)«  ea(x)+  b(x)I(\/2bU))  (4) 
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Ckn/  \  [M&O/N^w/JtSU/^U)] 
SNR(“^  =  1  +  SM/N^u)  +  S(aJ)/N2(uj 


(7) 

(8) 


We  have  also  included  Appendix  A  to  demonstrate  that  the  expression  on 
the  right  hand  side  of  Eq.  (4)  is,  in  fact,  a  lower  bound  on  Pe(x).  Hence,  by 
substituting  Eq.  (4)  into  Eq.  (3)  the  inequality  sign  is  preserved. 

It  is  important  to  observe  that  since  Eq.  (7)  uses  arbitrary  spectral 
functions,  the  lower  bound  can  be  applied  to  investigate  a  wide  class  of 
signals.  As  pointed  out  before,  the  study  will  be  separated  into  two  parts: 
in  Section  II  we  consider  baseband  signals.  Analytically  it  appears  to  be  a 
simpler  case  and  should  therefore  be  understood  first.  In  Section  III  we 
consider  bandpass  signals.  In  that  context  all  the  results  concerning 
narrowband  signals  [9]  will  be  included  and  referred  to  as  a  special  case. 


II.  BASEBAND  SYSTEMS 


In  this  section  we  shall  concentrate  on  signals  whose  power  is 
distributed  about  zero  frequency.  To  simplify  the  form  of  results,  let  us 
consider  the  following  special  case 


S(«) 


(u>  [  <  W/2 


(9) 


0 


u|  >  W/2 


b(x) 


SNR  sin^u)x/2 
1  +  SNR  si n^u>x/2 


du 


(13) 


To  obtain  the  lower  bound  one  must  substitute  Eqs.  (12)  and  (13)  into 
Eq.  (4)  and  Eq.  (4)  into  Eq.  (3)  successively,  and  carry  out  the  indicated 
algebraic  operations.  Since  we  are  primarily  interested  in  the  ambiguity 
phenomenon  we  shall  assume  that  WD/2*  >>  1,  say,  by  at  least  a  factor  of  3. 

In  that  case,  following  some  rather  extensive  algebra  manipul ations  outlined 
in  Appendix  B,  it  is  shown  that  the  lower  bound  consists  of  essentially  three 
disjointed  segments  as  suggested  by  the  following  equation: 
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1D2/12  (WT/2tt)SNR  <  a 

THRESHOLD  a  <  (WT/2w)SNR  <  6  (14) 

24*/W3TSNR  (WT/2tt)SNR  >  6 

where  in  the  threshold  region  the  lower  bound  varies  essentially  exponentially 
with  (WT/2tt) .SNR.  This  result  is  illustrated  in  Figure  4.  The  lower  limit  of 
the  threshold  region  is  given  by 

a  =  0.92  =  -0.36  db  (15) 

The  upper  limit,  8,  is  the  solution  to  the  following  transcendental 
equation 

(b/2)J1b/2)  =  (6/WD)2  (16) 

Note  that  Eq.  (16)  has  two  solutions.  We  only  consider  the  larger  one. 
In  Figure  5,  8  is  plotted  and  tabulated  as  a  function  of  WD/2u  for  the 
convenience  of  the  reader.  Thus  in  a  logarithmic  (db)  scale  8  'aries  from 
about  13  db  for  moderate  WD  products  to  about  16  db  for  exceedingly  large  WD 
products. 

To  put  these  results  into  perspective,  we  first  observe  that  a  m.s.e.  of 
0/12  corresponds  to  a  random  variable  uniformly  distributed  in  [-D/2, 

0/2].  Such  a  performance  level  can  always  be  achieved,  regardless  of  signal 
observations.  The  first  line  in  Eq.  (14),  therefore,  reads  that  for 
combination  of  WT/2-rr  and  SNR  such  that  their  product  (the  so-called 
post-integration  SNR)  does  not  exceed  ~C  db,  signal  observations  are 
completely  dominated  by  noise  thus  essentially  useless  for  the  delay 
estimation. 

We  next  observe  that  m.s.e.  predictions  based  on  the  Cramer-Rao 


inequality  yield  the  following  lower  bound  ([9],  Appendix  B). 
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e 


> 


(17) 


Substituting  Eq.  (10)  into  Eq.  (17)  and  carrying  out  the  indicated 
integration  one  immediately  obtains  the  third  line  in  Eq.  (14).  Hence,  if  the 
post-integration  SNR  exceeds  s,  the  ambiguities  in  the  differential  delay 
observations  can  essentially  be  resolved.  Cross  correlating  the  received 
signals  one  obtains  a  differential  delay  estimate  whose  m.s.e.  is  closely 
characterized  by  the  CRL8. 

In  the  derivation  of  Eq.  (14)  we  have  assumed  that  WO / 2-rr  >>  1.  This 
condition  means  that  the  correlation  time  (inverse  bandwidth)  of  the  signal  is 
small  compared  with  the  maximum  expected  receiver-to-receiver  delay,  or  that 
the  spacing  between  receivers  is  large  compared  with  half-wavelength  of  the 
modulation  frequency  (signal  bandwidth).  Only  in  that  case  we  are  dealing 
with  the  possibility  of  a  significant  ambiguity  problem.  There  remains  the 
question:  How  the  results  stated  by  Eq.  (14)  are  affected  when  this  condition 
is  not  satisfied?  In  Figure  6  the  lower  bound  is  generated  numerically  (by 
exact  integration  of  Eq.  (3))  and  plotted  as  a  function  of  the 
post-integration  SNR  (in  a  logarithmic  (db)  scale)  for  fixed  WT/2it  and 
different  values  of  WD/2it.  The  threshold  points  a  and  a  were  calculated  using 
Eqs.  (15)  and  (16).  Note  that  is  normalized  by  D^/12  so  that  we 
actually  measure  the  relative  efficiency  of  the  attainable  m.s.e.  to  the 
a-priori  variance.  One  observes  a  significant  threshold  effect  whenever  WD/2ir 
»  1.  In  case  WD / 2 tt  <  1  (i.e.  in  ambiguity-free  situation),  one  observes  a 
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rather  smooth  transition  from  the  CRLB  to  the  D  / 12  asymptote.  Note,  in 
passing,  the  close  agreement  between  Eqs.  (15)  and  (16)  and  the  exact  limits 
of  the  threshold  region. 

We  finally  note  that  in  the  derivation  of  Eq.  (14)  it  is  assumed  that 

the  source  signal  and  the  additive  noises  are  spectrally  flat  over  the 

receiver  frequency  band.  Spectrally  flat  signal,  however,  has  a  highly 

oscillatory  correlation  function  indicating  a  serious  ambiguity  problem.  For 

signal  spectra  whose  correlation  function  (inverse  Fourier  transform)  is 

smoothly  varying,  the  ambiguity  phenomenon  may  not  be  as  critical.  The  lower 

2 

bound  in  that  case  is  still  composed  of  the  D  / 12  segment  and  the  CRLB,  but 
the  threshold  effect  is  expected  to  occur  at  a  lower  SNR.  This  effect  can  be 
studied  by  analyzing  the  lower  bound  for  various  signal  and  noise  spectra. 


III.  BANDPASS  SYSTEMS 


In  this  section  we  concentrate  on  signals  whose  power  is  distributed 
about  some  center  frequency  To  simplify  the  form  of  results  let  us 
consider,  in  complete  analogy  with  the  baseband  case,  the  following  special 
case 


S(u) 


!u  *  u0[  <  W/2 
ju  *  u0|  >  W/2 


(18) 


We  shall  further  assume  that  the  additive  noise  components  are 
spectrally  flat  over  the  signal  frequency  band  so  that  Eq.  (7)  assumes  the 
form 
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I  SNR  |u  *  u  |  £  W/2 

(19) 

0  lu)  ±  a)  |  > 'W/2 
1  o' 


where  SNR  is  defined  in  Eq.  (11).  Substitution  of  Eq.  (19)  into  Eqs.  (5)  and 
(6)  immediately  yields 


u0+W/2 


a(x)  = 


SNR  sin  wx/2)du 


VW/2 


uj0+W/2 

b(x)  =  ^  / 
VW/2 


SNRsin^tux/2  ^ 

- - - x -  U(j) 

1+SNRsin  wx/2 


(21) 


The  lower  bound  is  now  obtained  by  successive  substitutions  of  Eqs.  (20) 
and  (21)  into  Eq.  (4)  and  Eq.  (4)  into  Eq.  (3)  and  carrying  out  the  indicated 
algebra  operations. 

Let  us  first  assume  that  W/u>  <<  1.  Since  we  are  concerned  with  the 

o 

joint  effect  of  envelope  and  phase  ambiguities  on  the  attainable  m.s.e.  (see 
Figure  3),  we  shall  further  assume  that  WD/2n  >  1.  In  this  setting, 
therefore,  we  are  dealing  with  narrowband  signals  and  very  widely  separated 
receivers.  Following  some  rather  extensive  algebraic  manipulations  outlined 
in  Appencix  C,  it  is  shown  that  the  lower  bound  consists  of  essentially  two 

« 

distinct  threshold  effects  dividing  the  entire  SNR  domain  into  three 
disjointed  segments  as  suggested  by  the  following  equation 


-  12  - 


ID2/12 
THRESHOLD 

12tt/W3TSNR 

THRESHOLD 

tt/u2WTSNR 


(WT/2tt)SNR  <  Y 
y  <  (WT/2ir)SNR  <  6 

6  <  ( WT/2-rr )  SNR  <  u 

u  <  (WT/2ir)SNR  <  n 

(WT/2tr)SNR  >  n 


where  in  the  threshold  regions  the  lower  bound  varies  essentially 
exponentially  with  (WT/2*)SNR.  This  composite  result  is  illustrated  in  Figure 
7.  The  various  threshold  points  will  be  defined  shortly.  We  further  observe 
that  in  the  narrowband  case  the  CRLB  (obtained  by  substituting  Eq.  (19)  into 
Eq.  (17)  and  carrying  out  the  indicated  integration)  is  given,  to  an  excellent 
approximation,  by 


>  ir/uQWTSNR  (23) 

Thus  if  (WT/2ir)SNR  >  n,  the  lower  bound  coincides  with  the  CRLB.  This 

is  the  ambiguity  free  mode  of  operation. 

If  s  <  (WT/2ir)SNR  <  n,  the  lower  bound  exceeds  the  CRLB  by  a  factor  of 
o 

12(uq/W)  .  In  this  region  signal  observations  are  subject  to  unresolved 
phase  ambiguities,  however,  useful  estimate  of  the  differential  delay  can 
still  be  obtained  from  the  envelope  of  the  cross  correlation  function.  Note 
that  this  segment  of  the  lower  bound  coincides  with  m.s.e.  predictions  based 

i 

on  a  simplified  version  of  the  Barankin  lower  bound  [11].  We  further  note 
that  the  3rd  line  in  Eq.  (14)  exhibits  exactly  3  db  (a  factor  of  2)  loss  in 
m.s.e.  relative  to  the  3rd  line  in  Eq.  (22).  This  result  indicates  the 
following:  the  3rd  line  in  Eq.  (14)  is  the  CRLB  based  on  the  received  signal 
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envelopes.  An  estimator  which  realizes  that  m.s.e.  consists  of  shifting  the 
received  signals  to  baseband  (perhaps  because  of  unresolved  phase  ambiguities) 
prior  to  the  cross  correlation  operation.  The  latter  scheme,  therefore,  is 
always  inferior  to  a  scheme  which  first  cross  correlates  the  received  signals 
and  then  uses  the  envelope  of  the  cross  correlation  function  for  the  delay 
estimation. 

Finally,  if  ( WT / 2ir ) SNR  <  y,  the  lower  bound  is  essentially  characterized 

2 

by  the  constant  level  of  D1 112.  This  is  the  noise  dominated  region  where 
signal  observations  are  subject  to  envelope  ambiguities  as  well,  thus 
essentially  useless  for  the  delay  estimation. 

The  lower  and  upper  limits  of  the  more  critical  threshold  are  given, 
respectively,  by 

v  =  (2.76/ir2)U0/W)2  (24) 

n  =  (6/ir2)  (uQ/W  )2[£-1(  W2/24o>2  )  ]2  (25) 

where  jT*(  )  denotes  the  inverse  of_J(  ).  Eqs.  (24)  and  (25)  depends  only 
on  wQ/W,  the  ratio  of  center  frequency  to  signal  bandwidth.  Eq.  (25)  is  of 
particular  interest  since  it  represents  the  minimum  amount  of  post-integration 
SNR  required  to  achieve  the  CRLB.  Thus,  for  example,  if  (dq/W  =  10  (i.e.  10 
percent  signal  bandwidth),  y  =  14.5  db  and  n  =  28.3  db.  If  uQ/W  =  100 
(i.e.  1  percent  signal  bandwidth),  y  =  34.5  db  and  n  =  50.8  db.  One  further 
observes  that  the  threshold  region  is  not  infinitely  small  as  may  be 
interpreted  from  the  analysis  based  on  the  Barankin  lower  bound  [11].  For  10 
percent  signal  bandwidth  it  is  a  segment  of  13.8  db.  For  1  percent  signal 


bandwidth  it  is  a  segment  of  16.3  db. 


The  lower  and  upper  limits  of  the  secondary  threshold  are  given, 
respectively,  by 


r  =  a/2  =  0.46  =  -3.36  db 

(26) 

6  -  6/2 

(27) 

where  8  is  defined  by  £q.  (16).  Hence,  the  secondary  threshold  is  shifted  by 
exactly  3  db  (a  factor  of  2)  relative  to  the  threshold  phenomenon  found  in  the 
baseband  case. 

The  analytical  result  illustrated  in  Figure  7  is  derived  under  the 
assumption  that  W/(dq  «  1.  In  that  case,  the  lower  bound  exhibits  two 
distinct  threshold  effects.  The  phase  ambiguities  and  envelope  ambiguities 
occur  at  essentially  disjointed  segments  of  the  SNR  domain  and  therefore,  in  a 
sense,  strictly  additive.  This  conclusion,  however,  applies  only  to 
narrowband  signals.  As  the  signal  bandwidth  increases  the  two  threshold 
regions  come  close  ( i . e. ,  6  converges  to  y)  and  the  two  ambiguity  phenomena 
are,  in  fact,  indistinguishable.  This  effect  is  illustrated  in  Figures  8-10 
where  the  lower  bound  is  generated  by  exact  numerical  integration  of  Eq. 

(3).  Note,  in  passing,  the  close  agreement  between  the  various  threshold 
points,  calculated  using  Eqs.  (24)  to  (27),  and  the  exact  limits  of  the 
threshold  regions. 


* 
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Narrowband  Systems 

It  is  not  difficult  to  relate  the  composite  result  given  by  Eq.  (22)  to 
the  analysis  of  narrowband  systems  carried  out  in  [9].  In  [9]  the  lower  bound 
is  derived  under  the  assumption  that  the  signal  bandwidth  is  so  narrow  that 
its  correlation  time  (inverse  bandwidth)  exceeds  the  maximum  expected  delay 
(i.e.  WD/2ir  <  1).  The  envelope-type  ambiguities,  therefore,  are  completely 
eliminated  from  considerations.  Instead  of  having  a  threshold  effect,  there 
is  a  smooth  transition  from  the  first  line  of  Eq.  (22)  to  its  third  line  as 
illustrated  in  Figure  11.  The  point  at  which  the  transition  occurs  (denoted 
by  p  in  the  figure)  is  obtained  by  simply  intersecting  the  two  indicated 
lines.  One  finds 

P  =  (18/tr2)/(WD/2tr)2  (28) 

With  this  modification  Eq.  (22)  reduces  to  the  result  derived  in  [9]. 

IV.  CONCLUSIONS 

This  is  the  second  part  of  a  study  which  deals  with  the  problem  of 
passive  time  delay  estimation.  The  focus  here  is  on  systems  employing 
wideband  signals  and/or  arrays  of  very  widely  separated  receivers.  A 
modified  (improved)  version  of  the  ZZLB  is  used  to  analyze  the  effect  of 
additive  noise  and  signal  ambiguities  on  the  attainable  mean  square 

i 

estimation  errors. 

We  first  concentrated  on  baseband  (lowband)  signals.  When  the  lower 
bound  is  plotted  as  a  function  of  SNR,  one  observes  a  distinct  threshold 
phenomenon.  Above  a  critical  SNR  the  lower  bound  approaches  the  CRLB.  This 
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is  the  small  error  region.  Below  the  threshold,  the  lower  bound  quickly 
approaches  a  constant  level  depending  only  on  the  variance  of  the  a-priori 
search  domain  of  the  delay  parameter.  In  this  region  signal  observations  are 
completely  dominated  by  noise,  thus  essentially  useless  for  the  delay 
estimation.  Information  concerning  the  threshold  points  and  the  various 
segments  of  the  lower  bound  is  given  by  Eq.  (14). 

Delay  estimation  using  bandpass  signals  is  analytically  a  more 
complicated  problem.  Here  the  lower  bound  exhibits  two  distinct  threshold 
effects  dividing  the  SNR  domain  into  three  disjointed  segments:  at  high  SNR 
the  lower  bound  coincides  with  the  CRLB.  This  is  the  ambiguity-free  mode  of 
operation  where  differential  delay  estimation  is  subject  only  to  local 
errors.  At  moderate  SNR  (between  the  two  thresholds)  the  lower  bound  exceeds 
the  CRLB  by  a  factor  of  12 W )  where  u>q  and  W  are,  respectively,  the 
center  frequency  and  signal  bandwidth.  In  this  region  the  ambiguities  in  the 
received  signal  phases  cannot  be  resolved,  however,  a  useful  estimate  of  the 
differential  delay  can  still  be  obtained  from  the  received  signal  envelopes. 
At  low  SNR  the  lower  bound  approaches  a  constant  level  depending  only  on  the 
variance  of  the  a-priori  search  domain  of  the  delay  parameter.  In  this 
region,  signal  observations  are  subject  to  envelope  ambiguities  as  well,  thus 
essentially  useless  for  the  delay  estimation.  Information  concerning  the 
various  threshold  points  and  the  various  segments  of  the  lower  bound  is  given 
by  Eq.  (22). 
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APPENDIX  A 


Lower  Bound* on 


The  binary  decision  problem  under  consideration  here  is 
given  by 


Ho  -  A  -  '  &T»  *  Ck 

H.  :  LZ  -  CT,  z  a4* 


The  log  likelihood  ratio  test  (log  LRT)  between  Ho 
is  defined  by 

l ,  Jc g  i<i!M 

C  PcHm) 


and 


H. 


(A- 2) 


where  £.(  1_  ]  H-i  \  is  the  conditional  probability  density  of  'T 
under  Hs  hypothesis,  and  £"  represents  the  data  vector  (e.g.  time 
samples  or  Fourier  coefficients  of  the  received  signals).  Assuming 
that  Ho  and  u>  are  equally  likely  to  occur  (i.e. 
p  P  hO  -  A,  ^  a  decision  rule  which  minimizes  the 

probability  of  error  compares  the  log  LRT  to  zero  threshold. 

If  c  PC  we  decide  on  Hi  ,  if  ^  0  we  decide  on  Ho 

Hence,  the  minimum  attainable  probability  of  error  is  given  by 


oO 

!  p  '/>.)<«/ 

) 

o 


4. 


I 

X 


J.  |  H 1 )  c(  J 


-oo 


(A-3) 


Let  ^5)  denotes  the  characteristic  function  associated 
with  .  p.  ancj  are,  by  definition,  a 
Fourier  transform  pair  satisfying 
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c*  +  'i+* 


I  H*)  = 


%a)e,  Ji 


(A-4) 


<-0*3°° 


y 


and 


C  + 


w  oO 


n 't I H.)  -  jjTJ"  J  <Ji 


-iJ 


(A-5) 


‘■T 


The  integrals  in  Eqs.  (A-4)  and  (A-5)  must  be  evaluated  at 
C0-o  and  C,  -  0  ,  respectively.  However,  for  future  manipulations, 

let  Co  be  a  real-valued  parameter  in  the  open  interval  (0,1). 
Similarly,  let  C,  be  a  real-valued  parameter  in  (-1,0).  In  that 


case,  substituting  Eqs.  (A-4)  and  (A-5)  into  Eq.  (A-3) ,  one  obtains 

C'*r  m 


^  a,  cw 


-- 

^  j 


£%< 6  2 


i  at  w 


^/S)7s|£  c 


■  r 

> 


,oO 


‘•T 


c,-r 


— o 


=  z  jr-i  [  ~  j  r.'j'.'ii/iJ  ds 

‘  .  -t  0^1  £  ,  -  'VaO 

0 


73-6  ) 


-r 

Note,  in  passing,  that  the  transition  from  the  first  version 
of  Eq.  (A-6)  to  its  second  version  can  be  carried  out  if  and  only 
if  Co1>0  and  C,  <0  .In  this  setting  Q  **  *  )  is 

expressed  in  terms  of  -4  *  Cy  {  •  Now,  by  definition 


yi(S)=  E  {  i  A  =  0,  4 


(A-7) 


where  £  ^  ^  denotes  the  statistical  expectation  of  the  bracketed 
quantity.  Substituting  Eq.  (A-2)  into  Eq.  (A-7)  and  observing  that 
£  r  /(-T)  (so  that  the  statistical  expectation  operation  can  be 
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performed  with  respect  to  the  probability  density  of  V*  )  >  one 

immediately  obtains 


(A-8) 


_-S 

[  I  H/»)  c(  f 


(A-9) 


We  shall  now  generate  the  data  vector  V~  by  Fourier  analyzing 
r*(t)  in  Eq .  (1).  Since  signal  and  noise  are  assumed  to  be  zero 
mean  Gaussian  processes  and  the  components  of  Y~  are  generated  by 
linear  operations  on  these'  time  functions,  has  a  multivariate 

Gaussian  distribution 


f-'H  I H ; ) - 

where 


i 


(A-10) 


v<<  r  )  tol  -r  -  0  4.  CA-U) 

^  >j  ^ 

and  £"*  denotes  the  conjugate  transpose  of  .To  obtain  < S  ) 

one  must  substitute  Eq.  (A-10)  into  Eq.  (A-8)  and  carry  out  the  indicated 
algebra  operations.  For  observation  time  T  large  compared  with  the 
correlation  time  (inverse  bandwidth)  of  signal  and  noise  (i.e.  vvT/^rr  J  ) 
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the  Fourier  coefficients  (at  each  receiver  output)  associated  with  different 
frequences  are  statistically  uncorrelated.  In  that  case  the  required 
computations  become  relatively  easy.  Details  are  contained  in  (  [9  ] 

Appendix  A) .  The  result  is 

:  rr  (a— 12) 

K 

where 


*  -  ^rrv</T 


(A— 13) 


and 


P  (lOvc,  x ).  =  S^R  (  lOv<) 


(A-14) 


The  index  X  varies  over  all  components  in  the  signal  frequency 
band.  Following  very  similar  considerations  for  one  obtains 

x) 

vs  <- 


^,'5)-  JT  +  X)j 


(A— 15) 


Note  that  ^,(5)  and  ^,(5)  depend  only  on  X  so  that 
O  (a^-'x)  -  P*  (X  )  .One  must  now  substitute  Eqs.  (A-12)  and 

(A-I5)  into  Eq.  (A-6)  and  carry  out  the  indicated  integration.  It  can  be 
easily  shown  that  all  the  poles  (singular  points)  of  \^/5fS)are  located 
on  the  real  axis  outside  the  interval  [0,1]  ,  all  the  poles  of 

are  located  on  the  real  axis  outside  the  interval  [-1,0]  .  Hence,  the 

first  integral  in  Eq.(A-6)  can  be  evaluated  for  any  value  of  in  the 

open  interval  (0,1)  without  affecting  the  desired  result.  Similarly,  the 
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second  integral  in  Eq.  (A-6)  can  be  evaluated  for  any  v&lue  of  Ci 
in  (-1,0).  We  shall  find  it  most  convenient  to  choose  Co-  V2  } 

We  shall  further  make  the  change  of  variables  5  s  \  ■*  ^  in  the 

integral  and  £  -  in  the  second  integral.  In  this 


first 


setting  one  obtains 


.  1  _L 

-2. 


j 

-oO 


-  J  '  rr  T/  4  /  J  w^mP-LO*  \)  ? 

^  (  '  '  -  ; 


n  i 


-oO 


1  ^ 
^  K 


7— 4“i‘ 

_  X  -  '  -  J  J  -  A  _  'J’  U-O  . 


c 


— c — 


P>  -H? ) 


This  is  an  exact  expression  for  ^  *  )  .  To  generate  a  lower 

bound  on  (P  we  shall  use  the  inequality  ^  Q  whenever 

y  ^  0  .  Thus 


„..v 


- 1 


^ r'LJx,*)  ^  '  Wpr  HP(wa) 


(A- 17) 


Substituting  Eq.(A-17)  into  Eq.  (A-16)  one  obtains 

OO 


?  [p  ivo  v,] 'J 


I 


03 


.  „<*>,£  L-*,tf‘*w_±L 


-  a 


U‘ 


-\±A- 

K  'I -M  (10^*1  J 


(/9  -1$) 
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where  we  defined 


CL(>  )  -  -  2  I'n  -  r[  a  )] 


(A-19) 


b  ('x  )  X 


r  (to^  x) 


(A-20) 


k  f  j  r(u jKjx) 

The  integral  in  Eq .  ( A— 18)  can  be  evaluated  analytically  using 

[12],  p.  338,  formula  3,466-1.  The  result  is  given  by 


CV*/  -»  ,  /  , - — 1  v 

x)  2  e  <p(\hb(y)  ) 


( A-2 1 ) 


which  is  the  desired  result.  Substituting  Eq .  (A-14)  into  Eqs.  (A-19) 

and  (A-20),  one  obtains 


a  >J  =  -  *»■  'wOv<)-S(  >->Z 


K  -J 

T"  -SA  fl  ftAl  m  )  •  S(/vT  tO  k  X/i. 

C  X^  :  Z.  - - - ; - - - 

V'  ^  ,S/V  d  W/,  x/z 


(A-22) 


(A-23) 


ATT 


4  'f-* 


&CO 


where  we  define  * XV/t  •  For  large  WT  product  and  smoothly  varying 
signal  and  noise  spectra  the  function  *SV/T(oo)  S(>'1tOX/i,changes  only 
insignificantly  over  the  frequency  increment  of  ,  so  that  the  sums  in 
Eqs.  (A-22)  and  (A-23)  can  be  converted  into  the  integrals  in  Eqs.  (5) 
and  (6),  respectively. 


Eq .  (A-21)  can  further  be  used  to  generate  weaker  lower  bounds 
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! 

on  Pe(x) ,  which  will  become  useful  in  the  proceeding  analysis.  We  shall 
start  from  Eqs.  (A-19)  and  (A-20).  Since 


P(ujw,  Q 

i  -» 


4 


Hu,.,,.) 


(A-24) 


and  since 


r(LO«,x)]- 


P  (t-0  *  j  X  ) 

}  A  PiUtjY) 


£  3 .  r  x) 


(A-25) 


It  immediately  follows  that 


a b'x)  i  i  x  *> 

* 


(A-26) 


ar.d 

co )  a  \>(y)  ^  -  i  r  *) 

( A- 2 7 ) 


where  2b(x)  and  [a(x)+b(x)]  are  the  terms  appearing  in  Eq.  (A-21). 
Replacing  the  first  term  by  its  upper  bound,  and  the  second  term  by  its 
lower  bound,  ^.(^Ocan  further  be  bounded  by 


PeM 


£ 


4>(i/c(>f  ) 


(A-28) 


where  c(x)  and  d(x)  are  given  by 
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C<>)  =  c3>2  P(cjkj  x)  -  rf  '5.£a/?(uj«)  $U,l(ui«x/)Jau3 

©O  * 

- 

Wr/ifi»  1  ^  St^%o  X/z)  cluj  (A-29) 

Jo 

d(*)  :  }  I  p  (t-O^X)  s  X  [_5/*rt(u}*)-5c^%J*>/aJj  c,u> 

K 

,oO 

37  — * f  £r  I  [WKM'S^wXAycLo  (A_30) 

■b 

Using  the  inequality  5 l  vi^ox/;  f  Ox)  and  d/X)  can 

further  be  bounded  by 


Px  )  ^  c1/1 

(A-31) 

^(x)  ^  d***4 

(A-32) 

where 

OO 

r 

t-  L 

C  '  CO  ,S  /V  ft.  ( )  cAL"C 

0 

(A-33) 

and 

d“  =  j  [u^Afi(CM)]  JlO 

• 

(A-34) 

Substituting  Cfr)  and  difxj 

by  their  upper  bounds,  one 

immediately  obtains 


P<(*)  >  £  <p(cx) 


(A- 35) 


-  26  - 


We  note,  in  passing,  that  the  lower  bound  given  by  Eq.  (A-35) 
is  useful  only  for  small  x,  where  the  inequalities  in  Eqs.  (A-31)  and 
(A-32)  are  tight. 
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Appendix  Analysis  Of  the  Modified  ZZLB  For  Baseband  Signals 

The  modified  ZZLB  is  given  by  Eq.  (3),  rewritten  here  for 

reference 

r 

i7  5  J  X6-[(fc-x)  f^oo]  cJx  (B-i) 

0 

The  function  assumes  a  simple  form  if 

^x)  contains  several  well-defined  peaks.  This  is  illustrated  in  Figure 
12  for  a  typical  baseband  case.  Thus,  if  X^X^X^...  are  the  local 
maxima  of  Pfy)  ,  then  at  each  segment  J  ,  vjl  the  function 
(j  jjD- *)  Pc  fx)j  is  closely  bounded  by 

OpD-x)  ?ei*y  %  (£>  -x*»)  l^'xv.)  (b-2) 

Note  that  since  x)  ^X)_,  is  a  non-increasing  function  of 

x,  Eq.  ( B— 2 )  holds  for  arbitrary  set  of  X~’S»  We  further  note  that  the 
local  maxima  of  f3(x)  are,  in  fact,  the  local  maxima  (ambiguity  points)  of 
the  cross-correlation  function.  These  occur,  approximately,  at  (see  Figure 
2). 

X*  ~  Yi  -  0J  4j  Xj  "  '  ( B-3 ) 

Pt(Xr>)>  required  for  the  computation,  of  Eq .  (B-2),  is 

closely  approximated  using  the  expression  in  Eq.  (4).  Following  the 
derivation  in  Appendix  A,  that  expression  is  shown  to  be-*a  lower  bound. 
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Thus 


P*<M>  S'*''4  blx’)  VGSiZn 


(£-4) 


where  C\fv^)  and  0^**1  are  obtained  by  substituting  *  s  yn  into  Eqs, 
(12)  and  (13),  respectively.  One  obtains 

W/i 


ac*>)  z  ~  ^  X>s  [d  +  Sa/? •  b^L  oo  v^/i  ]  d<o 


vj  r  4 


WV./], 


WXk  /« 

4  +  i^R  $  ^S\  )c^  A 


and 


w/j. 


b(^v)  =  un 


S a/  d 

1  -1  St'*-1  u)  x*,/x. 


to 


( E-5 ) 


_  w  r  4 


^•TT  w>^/z 


/  w  X^/m 

«£/V  ^  ■  Svy«V-fX 
1  +  ^  -SOv^-A. 


d-a 


(B-6) 


where  ^ /T/j  .  Since  the  integrands  in  Eqs.  (B-5)  and  (B-6)  are 

even  and  periodic  functions  of  J1  with  a  period  of  TT  ,  and  since  we  are 
integrating  over  Fo^nTT/z^  *  these  equations  become  independent  of  n, 
i.e.  aftn)  ~Cx  and  b(V^  =-  b  ,  where 

JA 

-£yy(  4  ■*  S&R  so^Lj\  )c*-TV.  -  ~  llIilEzji  (e-7) 


a*  :  - 


Wf 


*rr  rr 


I  v/l 

WT  I  SMt.S'yy}-^ 


,  _ _  _  Ur  1 

b  '  J?f  tt  j  )  -»  ^  rr  ^  -*  S/vfi' 


(B-8) 
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The  integral  appearing  in  Eq.  (B-7)  can  be  found  in  [12], 
p.  594  formula  4.399.  The  integral  appearing  in  Eq.  (B-8)  can  be  modified 
to  a  form  found  in  [12],  p.  152,  formula  2.562.  Thus 


Pt  ('x*.')  > 


e*" b  ) 


(B-9) 


Substituting  Eqs.  ( B— 3 )  and  (B-9)  into  Eq .  (B-2),  one  obtains 

" *  j  ^  (D "b*  )  * »*  (B-10) 

Since  3  -X>,  ^  L-iT/w -X  for  all  X  €  [_Xx  “/Ih/O,  knj  ,  it  immediately 
follows  that 

(p  \£o  )  <•  IWTVv  (B_n) 

“  '  0  D-lTT/w  ^  ^  D 

The  first  line  in  Eq .  ( B— 1 1 )  is  illustrated  by  the  dashed  line 
in  Figure  12.  The  second  line  in  Eq.  ( B— 1 1 )  indicates  that  zero  is  a 
better  bound.  For  'WJ>/*TT  »  {  ,  the  lower  bound  presented  by  Eq.  (B-ll) 
is  very  tight  except  for  values  of  X  *n  the  vicinity  of  X  *  0,  where 
changes  rapidly. 

To  take  this  effect  into  account,  we  observe  that  for  small  x, 
P<(X)  is  closely  bounded  using  Eq  .  ( A— 3 5 ) .  In  that  region,  therefore, 

Pg(xJj  is  closely  bounded  by 


■» 


-4V 

-  (  !>-*)£  <j)(cx) 


(B-12) 


where  in  the  transition  from  the  first  version  of  Eq .  (B-12)  to  its  second 

version  we  observe  that  f*)J  *  j^x)  whenever  )  is  a  non-increasing 
function  of  x.  The  parameters  c  and  d  are  obtained  by  substituting  Eq . 

(10)  into  Eqs.  (A-33)  and  (A-34),  respectively.  One  obtains 


c  '  htt  ‘O  S**'  - ajTrr — 


(B-13) 


W/  1 

6^TT  |  <^0O  - 


VO*  r  SfiQ1 
S’.  A"  -TT 


(B-14) 


For  future  manipulations  we  note  that 


( B— 15) 


Combining  Eqs.  ( B— 11)  and  ( B-12)  ,  £-[CD-x)  P*.  )J  is  tightly 


bounded  by 


Cl  +b  _ ,  ^ 

Pft6< )]  ^  ■  (D-ai'/w-x)  €  ^  4y 


J5  -ZTT/w  <  *  £  s 


where  \  is  defined  by 


(B-17) 


Note  that  J  Is  approximately  the  point  at  which  the  first  and 
second  lines  in  Eq.  ( E— 1 6 )  intersect.  Substituting  Eq .  (E-16)  into  Eq  . 
(B-l)  one  obtains  the  following  lower  bound 


,  ^  _dv  r^/w  ***> 

Zijx/D-  <p(c  xjdix  /(r)-iTT/w-A)dxj€  (fi'JTb)  ( B-l  8) 

l 

Carrying  out  some  straight  forward  algebra  manipulations,  it  can 
easily  be  shown  that  ^.Tt /\~J  .  Since  we  are  assuming  that  1>  , 

the  first  term  on  the  right  hand  side  of  Eq .  ( B— 1 8 )  is  closely 


-D  -xv/yj 


approximated  by 

r  2 

5  j  x(  D~y)&  $  (c  x)c*x 


c 


xe  <fi'c>)dx 
JHy 


*  =  c‘j  y-^-^rB)l  ^  (IM9) 


where  in  the  transition  from  the  second  versions  of  Eq.  (B-19)  to  its 
third  version  we  have  substituted  Eqs.  (B-15)  and  (B-17).  Similarly,  the 
integral  appearing  in  the  second  term  of  Eq.  (B-18)  is  closely 
approximated  by 

J)-A7T/w 

^  J  ->)c(.x  &  lf/6  ( B-20) 

l 

Substituting  Eqs.  (B-19)  and  (B-20)  into  Eq .  (B-18),  one 


obtains 
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■&yp[- 


cMb 


-io  (wr/jjr) 


|  4  r  £  J 


( E-2 1 ) 


Numerical  integration  indicates  that  Eq.  (B-21)  is  an  excellent 
approximation  to  the  exact  result  (Eq.  (B-l)).  The  lower  bound  presented 
by  Eq.  (E-21)  exhibits  two  distinct  asymptotes:  As  ,  the 

variables  a,  b  and  c  approaches  zero.  In  that  limit  the  second  term  in  Eq . 
(B-21)  becomes  the  dominant  term  and  the  lower  bound  approaches 


?  r>Vu 


(E-22) 


As  ,  a+b  - b  -oo  and  b  — >  wt/mit  .  In  that  limit 

the  first  term  in  Eq.  (B-21)  becomes  the  dominant  term  in  the  sum  and  the 
lower  bound  approaches 


V  wvTrr1 


6a  ^ 


j. 
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Q  )  etc 


wt/jitt  »  1 


ClJ  > 
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( B-23) 


One  immediately  identifies  Eq .  (B-22)  with  the  first  line  of 

Eq.  (16),  and  Eq.  (B-23)  with  its  third  line. 

The  transition  from  the  £>l/\ 2  asymptote  to  the 
asymptote  essentially  starts  when 


( E-24) 
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and  is  essentially  completed  when 

-f  - 1/** 


where  in  the  transition  region  (the  so-called  threshold  region)  the  lower 
bound  varies  essentially  as  <p(vllb' )  *  One  must  now  substitute  Eqs. 

(B-7),  (B-8)  and  (B-13)  into  Eqs.  (B-24)  and  (B-25),  and  solve  these 
equations  with  respect  to  SNR  in  order  to  obtain  the  corresponding  3  db 
points.  Suppose,  for  the  moment,  that  these  solutions  are  obtained  at 
SNR  ^ <  1  .  In  that  case,  Eqs.  (B-7)  and  (B-8)  can  be  approximated, 

without  incurring  any  significant  errors,  by 


a  c 


WT  £/V/? 
j.rr  H 


( B-26) 


b  S' 


w t  a 

*rr'  a. 


( B-27 ) 


With  these  approximations ,  Eqs.  (B-24)  and  (B-25)  assume  the 
simplified  forms 

7k  )  ~  4/M  (  B-28) 


and 


R/A.  (fcfJtf/i  )  -  (6/^$  )l 


( B-29 ) 


where  R  is  the  post-integration  SNR  defined  by 


R  -  (vjT/^n)  S//fl 


( E-30) 


Denoting  by  #  =  o(  the  solution  to  Eq.  ( B— 28 )  and  substituting 
fl-lAO]1  -  »  one  immediately  obtains  Eq.  (15).  Denoting  by 

0,-fk  the  solution  to  Eq .  (B-29),  one  immediately  obtains  Eq .  (16).  We 


further  observe  that  in  terms  of  SNR  the  solution  to  Eq.  (B-28)  reads 
SNR  :=</(  w/T/ili)  ,  the  solution  to  Eq .  (B-29)  reads  SNR  -  (l  /(  • 

Thus,  for  WT/^i-n  »  {  ,  the  simplifying  approximations  made  in  Eqs. 


(B-26)  and  (B-27)  may  affect  these  results  only  insignif icantly . 
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Appendix  C .  Analysis  Of  The  Modified  ZZLB  For  Bandpass  Signals » 


The  modified  ZZLB  Is  given  by  Eq.  (3),  rewritten  here  for 

reference . 


1 

D 


o 


( C—  1 ) 


For  values  of  x  in  the  vicinity  of  vro  ,  the  function 
P^(*)  is  closely  bounded  using  Eq .  (B-12),  where  the  parameters  c 

and  d  are  obtained  by  substituting  Eq.  (19)  into  Eqs.  (A-33)  and  (A-34), 
respectively.  Assuming  that  >>  i  (so  that  one  can  ignore  terms  on 

the  order  of  'v/Vo.,)  relative  to  1),  C1  and  d  are  given,  respectively, 
by 


d1 


WT 

a  TT 


(C-2) 


Vv 


TT 


LOo  &  /*  Q. 


( C-3 ) 


For  future  manipulations  we  note  that 


(c/d)H  -  Hwr/^)  <C-4) 

For  values  of  x  away  from  x~0  ,  a  tighter  lower  bound  on 
C-[(D~>0  lx-]}  can  be  generated  from  the  local  maxima  of  ^(*1  .  This  is 

illustrated  in  Figure  12  for  a  typical  baseband  case.  In  the  bandpass 
case,  however,  PeG*'  has  two  sets  of  local  maxima.  One  set,  X»"(<An/w)^ 
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yy;  4tJL  ...  .  is  associated  with  the  ambiguities  in  the  envelope  of  the 

cross-correlation  function,  the  other  set,  =  Lj.  TyCo*)  />  ,is  associated 
with  the  ambiguities  in  the  phase  of  the  cross-correlation  function. 

We  shall  assume,  without  any  significant  loss  in  generality  that 
Ooa/vj  K  is  an  integer.  In  that  case,  using  the  6et  and 

following  the  same  considerations  outlined  in  Appendix  B,  &jp' *)  Pt(*l  i6 
bounded  using  Eq .  (B-ll),  where  a  and  b  are  given  by 


w  r  /  I  ^  /T-»  3  *  <?  ' 
rr  ^  I  ' 


( C-5 ) 


WJ  -  1  ( C-6 ) 

b  JT- * 

Note  that  there  is  a  factor  of  2  difference  between  Eq (C-5) 
and  Eq.  ( B-7 ) ,  and  between  Eq .  (C-6)  and  Eq .  (B-8). 

We  shall  now  use  the  set  X»  to  generate  another  lower 

bound  on  G-[^  £>-*)  (*■!_  using  similar  considerations.  Making  use  of  Eq . 

(A-28),  Pe(Xr.}  is  bounded  by 

p,(t)  >  4>(£i?J )  ‘c-7> 


where  £ly„)  and  are  obtained  by  substituting  Eq.  (19) 

(A-29)  and  (A-30),  and  calculating  these  functions  for  X  -  V*  . 

some  straight  forward  algebra  manipulations,  one  obtains 

A 

S/iQ-  C  wL.a  ci  -xv 


into  Eq . 
Following 


( C-8) 
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(C-9) 


Using  the  inequality  £  w*1  S*.  £  SX  }  C(x^)  and  ^(X^)  can  be 


bounded  by 


.  *d  2, ,  ^  L 

<(**)<-£  Y» 

C[  G(  J  £  Vr 


where 


.  VJ^T  ,->A 

C  '  rr 


w"  T 

T2^T~ 


For  future  manipulations  we  note  that 


(C-10) 


(C-ll) 


(C-12) 


(C-13) 


( c  /a  f  -  3-°  ( w-r/iir) 


(C-14) 


Suhst i t ut 1 ne 
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(C-15) 


Now  since 
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ptfr)]  2  (.©  '?»)  P4(V«) 


Vv  -aTV^o  ^  £  y^ 


(C-16) 


and  since  tve  right  hand  side  of  Eq.  (C-15)  is  a  monotonically  decreasing 
function,  it  immediately  follows  that 

Pefc)j>  (  D-^ir/w0- x)€  (X-JTT/Uo)]  (C-17) 


Combining  Eqs.  (B-12),  (C-17)  and  ( B— 1 1 )  (in  the  given  order), 


Cr[(J>-^  %(*)}  is  tightly  bounded  by 


!  -  3iu  *  jrr/u;,)* 

^  '  J  ^  «>18> 
i  &1  I  ^  VC  *  *  ^  -  ^■'r/kv 

:  0  3  -j?tvW  $  *  <-  D 


Note  that  is  approximately  the  point  at  which  the  first  and 
second  lines  in  Eq.  (C-lf)  intersect,  is  approximately  the  point  at 

which  the  second  and  third  lines  intersect.  Substituting  Eq .  (C-18)  into 

Eq.  (C-l),  one  obtains  the  following  lower  bound 
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3  j  >(£>-J.TT/w- /)^V 


e*iL4>(6V) 


( C-2 1 ) 


Following  some  straight  forward  algebra  manipulations,  it  can 
easily  be  shown  that  J,  •d^^n'/Lo,  ,  and  that  ^4JT7>W  .  Since  we  are 

assuming  that  < <  JO  ,  the  lower  bound  can  be  approximated,  without 

incurring  any  significant  errors,  by 

rl'  -cH/4* 
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where  in  the  transition  from  the  second  version  of  Eq.  ( C— 2 2 )  to  its  third 
version  we  have  substituted  Eqs.  (C-4),  ( C— 14 ) ,  (C-19)  and  (C-20). 

Substituting  Eqs.  (C-2),  (C-12),  (C-5)  and  (C-6)  into  Eq. 
(C-22),  one  can  now  generate  m.s.e.  predictions  for  any  pre-specif ied  SNR. 
We  further  note  that  the  second  and  third  terms  in  Eq.  (C-22)  contribute 
significantly  to  the  sum  only  when  SNR  4.4 f  .  In  that  region,  a  and  b  can 
be  approximated,  without  incurring  any  significant  errors,  by 


w  r 

TT  '  H 


(C-23) 


VvT 

xxt  X 


(C-24) 


With  these  considerations,  the  lower  bound  assumes  the  form 

..Tf-,'*' 

^ p IT 


( C-2 5) 


where  R  is  defined  as  the  product  of  VvT/jjr  by  SNR  (the  so-called 
post-integration  SNR). 

We  shall  now  examine  few  limiting  cases:  If  R  f  ,  the  third 
term  on  the  right  hand  side  of  Eq .  (C-25)  becomes  the  dominant  term  in  the 

sum  and  the  lower  bound  approaches 
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If  4  the  second  term  becomes  the  dominant  term  in 


the  sum  and  the  lower  bound  approaches 
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(027) 


Finally,  if  P  l/^  (vj,/w)  }  the  first  term  in  the  sum  becomes 

the  dominant  one,  and  the  lower  bound  approaches 
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xM 


WT/>7T»  1 


^3* 


1- 


-Uui  K 


(028) 


One  immediately  identifies  Eqs.  (C-26),  ( C— 27 )  and  (028)  with 
the  first,  third  and  fifth  lines  of  Eq.  (22)  respectively. 

The  transitions  from  the  first  line  of  Eq.  (22)  to  its  third 
line  essentially  starts  when 
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4)'vsr)  =  % 


( C-29) 


and  is  essentially  completed  when 


£  ^  >  -  d 


(wV6)R 


(C-30) 


where  in  the  transition  region  the  lower  bound  varies  essentially  as 
.  Denoting  by  X  the  solution  to  Eq .  (C-29)  and  using  the 
similarity  between  Eq .  (C-29)  and  Eq .  (B-28),  one  immediately  obtains  Eq . 

(26).  Denoting  by  the  solution  to  Eq .  (C-30)  and  using  the 

similarity  between  Eq.  (C-30)  and  Eq .  (B-29),  one  immediately  obtains  Eq 


(27). 


The  transition  from  the  third  line  of  Eq .  (22)  to  its  fifth 


line  essentially  starts  when 


—  *  K  j  r  Vm 


J  <9  L 


and  is  essentially  completed  when 


(C-31) 


w l/ il')  R 


( C-32) 


where  in  the  transition  region  the  lower  bound  varies  essentially  as 
(f)  ^  nH  )  •  Denoting  by  ({: and  ({  =•  j?  the  solutions  to  Eqs 

(C-31)  and  (C-32)  and  following  some  straight  forward  algebra 
manipulations,  one  immediately  obtains  Eq.  (24)  and  (25),  respectively. 


-  43  - 


References 


1.  V.H.  MacDonald  and  P.M.  Schultheiss,  "Optimum  Passive  Bearing 
Estimation  in  a  Spatially  Incoherent  Noise  Environment",  J.  Acoust. 
Soc.  Am.  46^,  pp.  37-43,  1969. 

2.  W.R.  Hahn,  "Optimum  Signal  Processing  for  Passive  Sonar  Range  and 
Bearing  Estimation",  J.  Acoust.  Soc.  Am.  _5£  (1),  pp.  201-20 7. 

3.  G.C.  Carter,  "Time  Delay  Estimation  for  Passive  Sonar  Signal 
Processing",  IEEE  Trans.  Acoust.  Speech,  Signal  Processing, 

Vol .  ASSP-29,  pp.  463-470,  June  1981. 

4.  Y.T.  Chan,  R.V.  Hattin  and  J.B.  Plant,  "The  Least  Squares  Estimation 
of  Time  Delay  and  Its  Use  in  Signal  Detection",  IEEE  Trans.  Acoust. 
Speech,  Signal  Processing,  Vol.  ASSP-26,  No.  3,  pp.  217-222,  1978. 

5.  W.R.  Hahn  and  S.A.  Tretter,  "Optimum  Processing  for  Delay-Vector 
Estimation  in  Passive  Signal  Arrays",  IEEE  Trans,  on  Information 
Theory,  Vol.  IT-19,  No.  5,  pp.  608-614,  1973. 

6.  B.V.  Hamon  and  E.J.  Hannan,  "Spectral  Estimation  of  Time  Delay  for 
Dispersive  and  Non-Dispersive  Systems",  J.  Appl.  Statis  2_3  (2), 
pp.  134-142,  1974. 

7.  C.H.  Knapp  and  G.C.  Carter,  "The  Generalized  Correlation  Method  for 
Estimation  of  Time  Delay",  IEEE  Trans.  Acoust.  Speech,  Signal 
Processing,  Vol.  ASSP-24,  No.  4,  pp.  320-327,  1976. 

8.  A.  Weiss  and  E.  Weinstein,  "Composite  Bound  on  the  Attainable 
Mean-Square  Error  in  Passive  Time-Delay  Estimation  from  Ambiguity  Prone 
Signal",  IEEE  Trans,  on  Information  Theory,  Vol.  IT-28,  No.  6, 

pp.  977-979,  November  1982. 

9.  A.J.  Weiss  and  E.  Weinstein,  "Fundamental  Limitations  in  Passive  lime 
Delay  Estimation  -  Part  I:  Narrowband  Systems",  IEEE  Trans.  Acoust. 
Speech,  Signal  Processing,  Vol.  ASSP-31,  pp.  472-486,  April  1983. 

10.  J.P.  Ianniello,  E.  Weinstein  and  A.  Weiss,  "Comparison  of  the 
Ziv-Zakai  Lower  Bound  on  Time  Delay  Estimation  with  Correlator 
Performance",  Proceedings  of  the  IEEE  International  Conference  on 
Acoustics,  Speech  and  Signal  Processing  (ICASSP  83),  April  1983. 

11.  S.K.  Chow  and  P.M.  Schultheiss,  "Delay  Estimation  using  Narrowband 
Processes”,  IEEE  Trans.  Acoust.  Speech,  Signal  Processing, 

Vol.  ASSP-29,  No.  3,  pp.  478-484,  June  1981. 

12.  I.S.  Gradshteyn,  I.M.  Ryzhik,  "Tables  of  Integrals,  Series  and 
Products",  Academic  Press,  1980. 


-  44  - 


Figure  Captions 

Fig.  1.  Typical  Narrowband  Signal  Cross  Correlation 
Fig.  2.  Typical  Easeband  Signal  Cross  Correlation 
Fig.  3.  Typical  Eandpass  Signal  Cross  Correlation 
Fig.  4.  Composite  Bound  on  <£■■*  -  Easeband  Systems 

Fig.  5.  ^  vs.  \?J ~ 

Fig.  6.  <£A  VS.  Post  Integration  SNR  for  Easeband  Signals. 
Fig.  7.  Composite  Bound  on  Eandpass  Systems 
Fig.  8.  Normalized  m.s.e.  vs.  Post  Integration  SNR  for 
Fig.  9.  Normalized  m.s.e.  vs.  Post  Integration  SNR  for 
Fig.  10.  Normalized  m.s.e.  vs.  Post  Integration  SNR  for 
Fig.  11.  Composite  Eound  on  <£•*  -  Narrowband  Systems. 

Fig.  12.  ^ 'V.)  and  6“£j.S>-y )  f°r  Typical  Easeband  Case. 
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